The 7 — )■ a phase transition in cerium is not isostructural 
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Using a small amount of radioactive ^^^Cd atoms implanted in cerium lattice to probe the hy- 
perfine electric quadrupole interaction, we have detected a quadrupole electron charge density com- 
ponent in a-Ce, which is absent in 7-Ce. The appearance of the quadrupole density in a-Ce pre- 
dicted by the quadrupole scenario of the phase transition demonstrates that the spacial symmetry 
is lowered at the 7 — >■ a phase transition. This finding makes cerium the first element where the 
symmetry change is driven exclusively by the valence electron degrees of freedom while the atomic 
centers of mass (cerium nuclei) occupy the face centered cubic positions. We also report on nuclear 
quadrupole interactions in other high pressure phases of cerium: a" (C2/m space symmetry) and 
a' (of-U structure). 

PACS numbers: 81.30.-t; 76.80.-|-y; 71.70.Jp; 71.27.+a; 64.70.K-; 61.50.Ks; 71.10.-w 



The crystal space group symmetry is the most impor- 
tant characteristic, which is used for the identification of 
a macroscopic collection of atoms or molecules. A tran- 
sition to another allotropic form of an element is there- 
fore fully described by the change of its space symmetry. 
This fact is exploited in the Landau theory of phase tran- 
sitions but the consideration fails at the apparently 
"isostructural" 7 a phase transition in cerium Q oc- 
curring at a pressure of 0.7 GPa at room ternperature 
and accompanied by a 16% volume collapse [1, To 
distinguish the essentials of these twin phases one has 
to introduce a very important characteristic equal in its 
significance to the crystal space symmetry. Since the es- 
tablishing of the face centred cubic (fee) symmetry of 
both phases in 1949 0] the transition has become a bat- 
tleground of various microscopic theoretical models offer- 
ing different parameters controlling the phase transition. 
Two models - a Kondo volume collapse [3, II] and a Mott- 
like transition for 4/ electrons - are well represented 
in the literature Despite much effort, there is still no 
consensus on the issue [9|. Needless to say, these models 
take the isostructural nature of the phase change as a 
well established fact. 

On the other hand, on the basis of thermal expansion 
and compressibility data it was suggested [l^l that in Ce 
instead of the critical point 0, @l 3' tricritical point |l| 
is realized. In such first order phase transition 

between two phases of different space symmetry contin- 
ues beyond the tricritical point as a second order phase 
transition 10]. Therefore, a-Ce should have a lower 
symmetry than 7-Ce [lo| . A mechanism of symmetry 
lowering was provided by a new scenario of quadrupo- 
lar electron charge density ordering llj, Il2| , which con- 
trary to the previous approaches predicts a hidden and 
very peculiar symmetry change at the 7 — >■ a transi- 



FIG. 1: Triple-(f antiferroquadrupolar (Pa3) structure of a- 
Ce proposed by Nikolaev and Michel ilV\. Quadrupoles rep- 
resent the I = 2 valence electron (4/ -|- 5id6s^) charge density 
distribution. 



tion. Fig. [TJ In the present work we present experimen- 
tal evidence for the presence of local quadrupolar elec- 
tronic densities in the a phase and the absence in the 7 
phase. The transition 7 ^ a is of the first order and 
is driven by the minimization of the electron repulsion 
between the neighboring cerium sites. The active elec- 
tronic mode belongs to the X point of the Brillouin zone 
and involves its three arms, Qx, (fx, that causes the 
quadrupolar order parameter components alternate sign 
along X, y, and z-axis. The space group symmetry low- 
ering, _Fm3m(7-Ce) — >■ Pa3 (a-Ce), is accompanied by a 
uniform lattice contraction so that the fee structure of the 
atomic centers of mass (cerium nuclei) is fully conserved 



11 



12|. The long range order is due to the orientation 
of the local 1^™^" quadrupolar charge density compo- 
nent of the valence electrons of cerium (4/ -f bd&s^) on 
four different sublattices. Fig. [TJ One can speak of a 
triple-(f antiferroquadrupole (3-q-AFQ) structure. The 
change can not be captured by usual x-ray diffraction 
experiments. Such a scenario is in agreement with previ- 
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ous experimental results, and in addition reconciles the 
7 — a transformation with the Landau theory of phase 
transitions. The significance of the latter fact should 
not be underestimated. In particular, it implies that 
apart from the specific symmetry change the 7 — a phase 
transition is not fundamentally different from the other 
cerium transformations: a — > a"(monoclinic C2/m sym- 
metry), a" — >■ a' (orthorhombic a-U space symmetry), 
etc. [1, It is also clear that the other theoretical 
models for the 7 — ^ a phase transition treat it as an ex- 
clusive case which is not immediately applicable to other 
phase transformations. Unlike these other approaches, 
the theory of quadrupole ordering lH 12 1 appealed to 
the experimental refinement of the electronic structure 
of a-Ce. 

The phase transition of the Fm3m Pa3 type is 
not uncommon in molecular crystals. It occurs in ful- 
lerite Ceo [l3l> nitrogen crystal (a-N2) [IB] etc. A hidden 
Fmim — >■ PnSm symmetry change was found experi- 
mentally in Np02 [16j. The Pnim symmetry 17| there 
is another triple-g structure, which differs from Pa3 in 
the way the threefold axes and electronic quadrupoles 
are distributed over four sublattices The problem 

with a-Ce is that the proposed space symmetry does not 
alter the fee structure of the cerium nuclei. Although 
in principle the change can be detected by synchrotron 
radiation experiments, in practice there are many seri- 
ous technical problems, which so far have obstructed the 
experiment. 

The existence of the quadrupolar order in a-Ce and 
its site symmetry can be tested by other experimental 
techniques. In particular, methods of nuclear physics 
are often used to describe the crystal electric field po- 
tential V and the local site symmetry by measuring 
the electric field gradient (EFG) at nuclear positions: 
Vij = d^V/didj = x,y,z). Being a symmet- 

ric traceless second-rank tensor, Vij can be diagonal- 
ized. Then it is fully determined by the largest prin- 
cipal component, Vzz^ and the asymmetry parameter 



,y = {V^^ ~ Vyy)/Vzz (0 < 77 < 1) [ig|. The electric field 
gradient is proportional to the nuclear quadupole fre- 
quency (QF) vq — eQVzz/h {Q is the nuclear quadrupole 
moment), which is customary used for characterization of 
its strength. Through the electric quadrupole hyperfine 
interaction the electric field gradient at a lattice site is 
directly experienced by a probe nucleus. Such nuclear 
quadrupole interactions in solids are exploited by many 
famous methods (for example, by nuclear quadrupole res- 
onance, Mossbauer spectroscopy, etc.). Below we use 
a less used techni que of time-differential perturbed an- 
gular correlations [li| (TDPAC) with the "iln/^^Cd 
probe atoms implanted in the cerium lattice It 
is worth noting that the TDPAC measurements can be 
performed with very few impurity nuclei and in contrast 
to the Mossbauer spectroscopy it is not restricted to low 
temperatures. In the best cases the method yields an ac- 



curacy approaching that of nuclear quadupole resonance. 

The method exploits an anisotropic intensity distri- 
bution of 7 rays emitted from an ensemble of the (in- 
termediate) excited nuclear states of probe ^"'^^Cd atoms 
(171-245 keV cascade). The cascade proceeds via the 245 
keV level with the half-life Ti/a = 84 ns, spin / = 5/2 
and the nuclear quadrupole moment Q — 0.83 b. The 
precession of nuclear spins in a static electric crystal 
field changes the alignment of the ensemble with time, 
which is observed by measuring the 7 ray anisotropy 
[see Eq. ^ below]. The initial alignment is prepared 
by populating the excited states with the electron cap- 
ture decay of the 2.8 d ^-'^^In isotope. The ^-'^^In activ- 
ity was produced via the ^^^Ag (a,2n) ^^^In reaction by 
irradiating a silver foil with the 32 MeV a-beam. Af- 
ter that the nuclear -'^^-'^In/-'^^-'^Cd probes were implanted 
in the cerium lattice by melting Ce powder (about 500 
mg) with a small piece of the irradiated silver foil ( <0.1 
mg) in a special chamber under pressure of 8 GPa^2o|. 
The TDPAC measurements were performed on polycrys- 
talline samples of cerium metal at room temperature us- 
ing a four-detector spectrometer equipped with a small 
size hydraulic four arm press with a capacity up to 300 
ton [2^ . The high pressure up to 8 GPa was produced in 
a calibrated "toroid"-type device with NaCl as pressure- 
transmitting medium, Ref. [llj. Non-hydrostaticity of 
the transmitting medium was checked by measuring the 
^^^Cd-TDPAC spectra of silver (fee) under high pressure 
and observed to be negligibly small. 

The angular time correlation is characterized by the 
second order perturbation factor of the static electric 
quadrupole interaction (QI) for the nuclear spin / = 5/2, 
which describes its precession due to a hyperfine interac- 
tion [1^: 

3 

G22(i; J^Q, ??, A) = S20 + XI cos(a;„t) e"'^""-*. (1) 

Here t is the time delay, 52,1 are the amplitude coeffi- 
cients, bJn are the angular precession frequencies {n = 
1, 2, 3 for / = 5/2) related to the energy differences be- 
tween hyperfine levels split by the QI. The frequencies are 
functions of the EFG Vzz or, alternatively, vq, and the 
asymmetry parameter 77. A is the half- width of the QF 
Lorcntzian distribution accounting for oscillatory damp- 
ing. 

Recording the delayed coincidence spectra at angles 
7r/2 and tt between detectors, N{7t/2, t) and iV(7r, t), one 
obtains the usual angular anisotropy 



, N{^, t) - N{^/2, t) 
'N{tt, t) + 2iV(7r/2, ty 



(2) 



which is the TDPAC spectrum. The pressure evolu- 
tion of the TDPAC spectrum of ^^^Cd in Ce is given in 
Fig. [5J Experiments have been carried out at room tem- 
perature and covered four phases of elemental cerium: 
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FIG. 2: Room temperature TDPAC spectra of "^Cd in cerium under pressure. (A) at normal pressure: 7-Ce, fee [Fm3m] 
structure, no quadrupole order; (B), (C), (D), (E) a-Ce, 3-q-AFQ structure at 1.8, 3.1, 3.9, and 4.4 GPa, respectively. (F), 
(G) a"-Ce [C2/m] at 5.3 and 6.2 GPa; (H) a'-Ce [a-U-structure] at 7.8 GPa. 



7, a, a" and a'. It can be shown that R{t) = 
— ^22<52G22(i), where G22{t) is the perturbation fac- 
tor, Eq. dU), Q2 ~0.80 is the solid-angle correction and 
A22 = —0.17 is the unperturbed angular correlation co- 
efBcient for the 7 — 7 cascade of ^^^Cd. The EFG param- 
eters are determined from a least-square fitting of the 
TDPAC spectra in accordance with Eq. ([Ij. 

The extensive TDPAC data are presently available for 
the ^^"'^Cd impurity, for which the EFG as a function of 
temperature and pressure has been determined in sev- 
eral rare earth (RE) metals [2^ ll^ ■ It should be noted 
that the majority of RE metals crystallizes in hexagonal 
structures implying a nonzero EFG. In particular, the 
double hexagonal close packed (dhcp) /3 phase of cerium 
has been thoroughly investigated by Forker et al. [1^ . 

In the cubic symmetry the first nontrivial electron 
charge density contribution at each site is described by 
the cubic harmonic K4 (i.e. I = 4) thus leaving no room 
for an EFG tensor: Vij — 0. Therefore, one should find 
I'Q ss and Vij « if both 7 and a-Ce are assumed to 
be isostructural. Instead, from our data for cubic phases. 
Figs. [2] and [H we find that the conclusion holds only for 
7-Ce, while in a-Ce vq ^ and Vij ^ signaling the 
appearance of a hidden quadrupole charge density com- 
ponent at the probe ^^^Cd atoms. 

We first describe the TDPAC spectrum of 7-Cc, 
Fig. IHA). At ambient pressure and room temperature 
the fitting gives vq =3.0(5) MHz (77 = 0). Notice, that 
our TDPAC spectrum of 7-Ce is very close to the one ob- 
served by Forker et al. (upper panel of Fig. 1 of Ref. [25l] ) 
with the same value of j/q . A small electric field gradient 



attributed to lattice defects and other impurities is found 
in all cubic metals and merely refiects the purity of the 
material. 

From the structural data we know that the TDPAC 
spectra at pressures 1.8, 3.1, 3.9, and 4.4 GPa should 
be identified as belonging to the a- phase, Figs.[2][B)-(E). 
They are consistent and from fitting we obtain = 
11(1) MHz (V^, =0.54(4) xl02iV-m-2) and a uniaxial 
local site symmetry (rj = 0), Fig. [3l Furthermore, the 
values of i^q and Vzz are close to the TDPAC parameters 
for /3-Ce (vq =12.5(7) MHz, V^z =0.61(4) x lO^^V-m-^ 
[2^), and other noncubic phases of cerium. Fig. [31 The 
observed values of QF and EFG for a-Ce unambigu- 
ously rule out the cubic symmetry and indicate a hid- 
den quadrupolar order in this phase. This experimental 
finding can be rationalized only by assuming a 3-q-AFQ 
order in a-Ce, Fig. [T] The only theory predicting such a 
symmetry chan ge f or cerium was developed by Nikolaev 
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and Michel 

The "^Cd TDPAC spectra at pressures 5 to 8 GPa are 
shown in Fig.[2jF)-(H). They show two step-wise changes 
characteristic of two other phase transformations. An 
increase in the quadrupole frequency at 5.3 GPa with- 
out the change of the asymmetry parameter [rj = 0) is 
induced by the modification of the Ce crystal structure 
from the simple cubic symmetry (a-Ce) to the monoclinic 
a" phase. A partial "softening" of the EFG and vq in 
a-Ce at 4.4 GPa is also an indication of a proximity to 
the phase boundary. According to the angle-dispersive 
powder diffraction technique [l3| the space symmetry of 
a"-Ce is C2/m. From Fig. [21 [31 one sees that the mon- 
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FIG. 3: Pressure dependence of EFG {Vzz, right scale) and 
the nuclear QI frequency {vq, left scale) of ^^^Cd implanted 
in cerium. *Data for /3-Ce are from Porker et al, Ref. [25i] . 



oclinic C2/m structure is preserved up to 7.5 GPa. At 
pressure 7.8 GPa the ^^^Cd TDPAC spectrum character- 
istics change again. Now, it is the orthorhombic (a-U) 
structure [13]. The assignment is corroborated by a large 
value of 77 = 0.52, typical for the a-U symmetry. (For el- 
emental uranium in the a-phase 77 = 1 [2^.) 

Fig. |3] presents the overall pressure dependence of QF 
and EFG in Cc. At the phase transition boundaries 
QF and EFG change discontinuously indicating the first- 
order character of the transitions. Earlier, the step-wise 
behavior of QF with changes of the crystal structure was 
observed for a number of rare earths 2^ 2^ ■ 

In conclusion, our TDPAC experiments (Fig. [3]) de- 
tect an appreciable EFG in a-Ce comparable with EFG 
for noncubic phases (/3, a"), which border a-Ce in the 
pressure-temperature phase diagram. This finding rules 
out the cubic symmetry in a-Ce and evidences in support 
of a 3-q-AFQ order suggested by Nikolaev and Michel 



There, the local site symmetry is trigonal (C3) with 



a quadrupolar electron charge density component 
oriented along one of the main cube diagonals. Fig. [T] 
With further pressure increase, the evolution of the TD- 
PAC spectrum follows other phase transitions in cerium: 
a — T' a" (C2/m) and a" — >■ a' (a-U structure). 
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